UK Patent Application n 9 ,GB ,id2 369 735 n 3 »A 


(43) Date of A Publication 05.06.2002 


(21) 

Application No 0029420.7 

(51) 

INT CL 7 


H03F1/32 

(22) 

Date ot riling oz.iz.zuuu 



(52) 

UK CL (Edition T ) 




H3W WULPR 

(71) 

Applicant(s) 


H3TT6E 


Roke Manor Research Limited 




(Incorporated in the United Kingdom) 

(56) 

Documents Cited 


Roke Manor, Old Salisbury Lane, ROMSEY, 


GB 22S4973 A WO 96/26624 A1 


Hampshire, SOS1 0ZN, United Kingdom 





(58) 

Field of Search 

(72) 

Inventor(s) 


UK CL (Edition S ) H3W WULPR 


Paul Demain 


INT CL 7 H03F1/32 


Ben Toner 


ONLINE:WPl,EPODOOJAPIO 


John Domokos 



(74) 

Agent and/or Address for Service 




Siemens Shared Services Limited 




Intellectual Property Department, Siemens House, 




Oidbury, BRACKNELL, Berkshire, RG12 8FZ, 




United Kingdom 




(54) Abstract Title 

A linearization means for a component (eg an amplifier) in which the inverse transfer function for a 
predistorter is adjusted in dependence on the l/P and O/P 

(57) A method of linearising the output signal of a component e.g. power amplifier having an input signal 
comprising; taking an initial estimated representation of the pre-distortion inverse curve, manipulating said 
representation; pre-distorting said input signal by applying values of said representation to the input signal; 
determining the error in terms of linearity between the input signal and the output signal; and repeating 
previous steps such that the error is minimised. 

Fig. 5 Digital PredistoriJon Block Diagram 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy 
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Fig. 1 
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Fig. 2 Feed Forward Linearisation Technique 
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Fig. 3 Adaptive Analog Predistortion 
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Fig. 4 AdaptiveDigital Predistortion 
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Fig. 6 
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Fig. 9 Special Contour Generation Example 
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METHOD OF LINEARISING A SIGNAL 

The invention relates to a method of linearising a non-linear device 
and has particular but not exclusive application to linearsed power 
amplifiers. The invention also relates to linearised power amplifiers 
embodying the method. It should be noted however, that the invention can 
be embodied in any system which has a non linear characteristic. 

Figure 1 shows a schematic graph illustrating the characteristic of a 
non linear system or component such as a power amplifier. The X-axis 
represents nominally the magnitude of the input to the power amplifier at a 
particular instant and resultant magnitude of the output is shown on the Y- 
axis. An ideal linear system is characterised by the complete line A. 
However, in practice, the real characteristics of an amplifier are likely to 
be non-linear as represented by the broken or dashed line B. As can be 
seen saturation occurs at high magnitudes of input such that increasing 
input does not result in a corresponding increase in the output signal. 

Known linearisation techniques applied to power amplifiers 
essentially consist of three parts: a non-linear device i.e. an RF amplifier; 
an output signal measurement subsystem; and a linearisation subsystem 
i.e. feed forward or predistortion. In order to appreciate the invention prior 
art linearisation methods are briefly described. 


One solution to the linearisation is by means of so called feed 
forward linearisation which is illustrated schematically in figure 2. Such an 
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implementation includes a cancellation loop, feed forward loop and a 
measurement subsystem; the latter being essentially a selective pilot 
receiver measuring the power spectral density (PSD) of the residual tone at 
the output. The key aspect of this method is that the PSD is measured and 
hence large dynamic range is achieved. The measurement subsystem is 
based on pilot signal and not on a graphical method. This solution 
achieves a very high linearity and excellent signal to noise ratios in excess 
of 70 dB. The disadvantage of the pilot is that there is poor correlation 
between the minima and linearity and ACPR (Adjacent Channel Power 
Ratio;. This can be improved by a spread pilot scheme but this is complex. 
Another disadvantage of feed forward technique is the poor efficiency due 
to coupler delay line losses and the need for an error amplifier. 

Another method of providing linearisation is that of adaptive 
analogue predistortion. Figure 3 shows a schematic representation of an 
adaptive analogue predistortion system. This shows power amplifier 1 the 
output of which is additionally input into a complex demodulator 
whereafter coefficients of a polynomial, representative of the inverse 
characteristic of the power amplifier, are computed. This polynomial is 
generated and applied to the signal which is input into the power amplifier. 
In this way the input signal is modified such that attenuated parts of the 
output signal are amplified and phase adjusted to compensate for the non- 
linear effects; the input signal to the power amplifier is predistorted. There 
are however, three problems associated with such predistortion. This 
system necessitates the use of a complex down-converter. The complex 
down sampler has less dynamic range than that of the PSD measurement 
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receiver because the full signal bandwidth has to be digitised 
instantaneously. This is an issue particularly for W-CDMA (Wideband 
Code Division Multiple Access) application which has higher bandwidth 
than an EDGE carrier. Also because the polynomial generator comprises 
analogue circuitry it is inherently imperfect especially for high order terms 
because of delay mismatch, linearity, etc. This is limiting the bandwidth 
and the quality of curve fit hence the linearity also. It has also been found 
that polynomial generation does not achieve good fit always and in some 
cases may become unstable and produce ripple effects through the various 
sections of the characteristics. 

Digital pre-distortion is also known. In prior art digital predistortion 
systems again a complex down-converter is required in the measurement 
subsystem. This restricts the dynamic range, especially for broadband 
systems. Figure 4 shows a schematic block diagram of a known digital 
predistortion system. In such prior art systems spline and rational 
polynomials are used for linearising the amplifier. The nature of pre- 
distortion demands a method of applying inverse transfer function 
characteristics to the amplifier source data. The pre-distortion has been 
applied in this method by means of determining the coefficients of a 
polynomial representing the inverse characteristic, and points on the curve 
are stored in Look-Up Tables (LUT) generator. However there are 
problems with achieving stable polynomials that correctly characterise a 
typical amplifier inverse transfer characteristic. 

There is not always a deterministic relationship between PSD and 
the characteristics of the non-linear device. Therefore the optimum pre- 
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distortion, which is the inverse of the non-linear characteristics, cannot 
always be readily calculated. 

It is an object of the invention to provide a method of accurately 
describing typical amplifier characteristics for any system whilst 
maintaining stability and practical implementation ability. It is further an 
object of the current invention to overcome the above mentioned problems 
and to provide a linearisation system which is fast and effective. 

The invention comprises a method of linearising the output signal of a 
component having an input signal comprising; 

a) taking an initial estimated representation of the component 
inverse transfer characteristics; 

b) manipulating said representation; 

c) pre-distorting said input signal by applying values of said 
representation to the input signal; 

d) determining the error in terms of linearity between the input 
signal and the output signal; and 

e) repeating steps b) to d) such that the error is minimised. 


Additionally the above method is useable with PSD measurement 
alone and advantageously does not require phase information. If the 
lineariser is capable of operating without the phase information of the 
signal, much larger bandwidth and larger signal to noise ratio can be 
achieved with the measurement receiver. 
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The optimum pre-distortion characteristics may be determined by a 
predefined search method. The iteration may be carried out graphically 
using NURBS (Non - Uniform Rational B-Splines). In most cases, a single 
error value is required for the curve fit, which represents the overall 
spectral re-growth due to the residual non-linearity. 

Preferably the area between the target power spectral density and the 
measured power spectral density is minimised across the full range of the 
spectrum. In a preferred embodiment the iteration is done so as to 
minimise the error at particular regions of the spectra. Such a weighting of 
the different regions of the spectra is also found to play a major role in the 
stability and the convergence of the loop. 

Other advantages of the invention and its various embodiments are : 
the error value to represent the whole spectra with appropriate weighting 
for best convergence and stability; the best possible noise immunity 
without compromising the spectral resolution is achieved, it is ensured that 
global minima is achieved without random search, and provide a single 
error value for the linearisation loop from PSD measurement of the output 
signal and ensure convergence to best Error Vector Magnitude 

The invention will now be described by way of example only and with 
reference to the additional following figures of which: 
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Figure 5 shows a block diagram representing an embodiment of the 
invention. 

Figure 6 shows a preferred method of determining the linearisation error 
between the input and output signals. 

Figure 7 show an example of how the inverse transfer function curve may 
be represented and manipulated. 

Figure 8a and b shows an alternative method of the representation of the 
inverse transfer function curve. 

Figure 9 shows a flow diagram showing a preferred method of 
determining the error. 

Figure 5 shows a schematic block diagram of the system according to one 
embodiment of the invention and comprises the Power Amplifier 1 . The 
signal from the power amplifier is fed into a Power Spectral Density 
measurement unit 2. This unit measures the power spectral density of the 
signal; figure 6 shows a representation of a target spectrum and a 
measured spectrum. Measured PSD, MPSD, is measured at the output of 
the power amplifier and includes pre-distortion and Target PSD, TPSD, is 
how the measured PSD should be once pre-distortion is optimized. 
The advantage of using power spectral density as a measure of the signal 
is that it is done in the frequency domain and is a relatively simple 
process. Data relating to PSD is input to an Error Calculator 3, which 
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calculates the difference between the measured spectrum and a target 
spectrum which has been pre-computed and stored in look-up format. 

The result is then input into a Graphic Control Unit 4 whose 
operation is at the heart of the invention. The Graphic Control Units 
function is to determine the fundamentals of the inverse transfer function 
(curve) which should be applied eventually to the input signal to minimise 
the distortion. The output of the GCU is to a curve fit and interpolation 
unit 5, which essentially computes and interpolates the basic 
characteristics of the distortion inverse curve, so as to provide a whole 
series of point values which are stored in the pre-distortion coefficient 
Look-Up Tables (LUT) 6 at the resolution of the signal to be pre- 
distorted. (The term "coefficient" here represents high resolution 
coefficients which are interpolated from the representation, which may be 
a series of just a few points, and is not to be confused with the term 
coefficients of a polynomial representing the inverse transfer function as 
in prior art systems). These pre-distortion coefficients are then fed back to 
a linearising processor which essentially uses the appropriate data from 
Look-Up Table to modify the input signal to the Power Amplifier. The 
LUT contains the pre-distortion coefficients for all possible input values 
at a given bit resolution and are addressed using the baseband input value 
signals. 

The essence of the invention is that the inverse transfer function i.e. 
the pre-distortion which is to be applied to the input signal is^represented 
as look up tables as in the prior art, but that rather than it being provided 
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by calculating a series of coefficients of a polynomial, it is represented by 
a series of control points, stored in buffer 7. (alternatively described as 
"handles")* i.e. a curve through the control points or generated from 
NURBS (see later) which can be manipulated. 

The operation of the graphic control unit and manipulation of the control 
points will now be described in more detail. The graphical control unit 
begins to move each control point of the complex transfer characteristics, 
following a predefined methodology, and feedback from the error 
calculator enables the graphic control unit to decide whether the control 
point or handle has been moved in the right direction. 

Example 1 

Figure 7 shows an initial set up of control points, A to E representing the 
inverse transfer characteristics (or inverse curve) In this example all the 
control points are set to 1, and thus the "curve" or characteristic is a 
straight line. When these characteristics are stored as pre-distortion 
coefficients, the curve fit and interpolation unit would store all the values 
in the look up tables as L When applied to the input signal it would leave 
the input signal unchanged. 

In next stage the last of the control points F is slid vertically up the y axis 
by a fixed amount to new position F' . The coefficients representing a 
curve passing through the control points i.e. representing the complex pre- 
distortion function are then recalculated and fed into the look up tables. 
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The linear processor applies these to the input signal to the power 
amplifier and then the error calculation is once again made. Should the 
error increase then generally speaking the control point has been moved in 
the wrong direction. The control points are then moved (or "slid") in turn 
in a predetermined order and between each movement the inverse transfer 
junction (pre~distortion curve) characteristics are re-interpolated to the 
resolution stored in the look-up tables. Again when applied to the power 
amplifier the feedback determines the error and a better position for the 
handle is calculated. This process is therefore iterative and eventually the 
optimum position of the control points representing a distortion curve is 
determined so as to minimise the error. 

The skilled person would realise that there are an infinite number of ways 
in which the control points may be manipulated. They could be 
manipulated sequentially in ascending or descending order or indeed more 
than one control point may me moved at the same time. Various 
mathematical techniques may additionally be employed to ensure that the 
optimum position of control points i.e. which is the optimum pre- 
distortion curve. 

The term making an "initial estimated representation", as in the main claim 
includes 

a) making some estimates of how much error is on the signal and 
subsequently fit an inverse compression curve to this data, and 
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b) simply initialising the control point buffers and coefficient LUTs to all a 
set value e.g. ones (i.e. no predistortion is applied) 

However the inventors have determined that the first method of estimating 
the initial distortion levels is prone to bias and the accuracy of the estimate 
is not particularly consistent. The other option of starting the algorithm 
with no predistortion produces a more consistent result, and furthermore it 
is not necessarily any slower to converge the PSD error function to a 
minimum 

Example 2 

The inventors have additionally determined that the invention can be 
further improved by representation of the inverse curve, i.e. the inverse 

UtansfenSRiB^gWS^'^rt series of basis functions; which can be 
advantageously manipulated before summing to produce a NURBS curve 
are particularly effective at implementing the above mentioned 

. embodiment. Figure 8 a shows a series of basis functions and figure 8b 
shows the interpolated transfer function they represent. The transfer 
function simply comprises the superposition of a number, seven in this 
example, of individual basis functions, A to G. 

p Each individual basis function can be manipulated in turn by 
scaling, and after each scaling, the error between the target and measured 
power spectral densities can be measured. Effectively the individual basis 
functions are weighted. 


In view of the substantial advantages of NURBS modeling seen, 
this technique's flexibility enables its use in other linearisation contexts. 
The number of weights correspond to the number of basis functions and 
the number of basis functions correspond to the order of the curve. The 
basis functions are generated recursively to fit any arbitrary curve 
description and are evaluated at the desired interpolation resolution. Once 
the basis functions have been calculated, real time shaping of the curve to 
fit any other arbitrary curve is performed by simply adjusting the weights 
of each basis function and summing them in a predefined fashion. The 
nature of the simplistic and fast adjustment of curve shape makes this 
method of curve generation ideal for representing amplifier inverse 
transfer functions in digital hardware. 

The advantages of a NURBS curve compared to a polynomial curve 
are the following: 

a) independently definable regions of the curve; the summation of 
the basis functions at any point on the x-axis will result in the 
corresponding y-axis value. As the basis functions only have a finite range 
within the overall curve range, regions of the curve can be adjusted 
independently without undesirable effects on other regions. 

b) the stability of the system is inherent, due to the independently 
definable regions of the curve. When an adjustment needs to be made to 
one region of the curve (e.g. the compression is to be hardened) it can be 
assured that other regions of the curve will not be affected by ripples 
caused by the change. This, of course, could not be guaranteed with a 
polynomial as certain curve shapes may fall out of the boundary of that 
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particular order of polynomial, thus creating a ripple effect throughout the 
curve. 

c) the accuracy to fit to amplifier characteristics; the variation of the 
weight of one basis function will only affect the overall curve within the 
region of that particular basis function thus providing the ability for a 
curve to show predominately 7 th order characteristics in one half of the 
curve whilst remaining linear in the other half, which would not be 
achievable with a single polynomial. This effect is particularly useful when 
characterising an amplifier with a very good linearity at low signal values 
and high order compression at larger signal values. 

d) reduced real time computation complexity; since the computation 
of the curve points is carried out by summing each weighted pre-computed 
basis function at a particular x value, the basis functions for a particular 
system can be recursively calculated offline and stored in look-up tables. 
Real time computation then only requires each basis function to be read at 
a particular point, weighted, and summed. 

In a preferred embodiment of the invention the inventors have determined 
a preferred way of generating the error which is to be minimised; it is a 
graphical method for providing error feedback for adaptive linearisation. 
The error is calculated by comparing the measured PSD point by point 
with a pre-calculated desirable spectral contours. The error is the total 
area, by which the contours are exceeded, therefore, 
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/ 

E = ^D(0 - r arg *r(/) 

PSD and Target are the normalized spectra points in dB and 7 is the 
length of the spectra vector given that the error is integrated over the entire 
span of the PSD, 1, maximum noise immunity is achieved. A set of 
spectral contours is generated representing the desirable profiles through 
the linearisation process. The first contour in the set corresponds to the 
spectra resulting from the original non-linear characteristics of the typical 
device. The last contour in the set, the target contour, corresponds to the 
ideal (target) spectra. The intermediate contours in the set provide a 
desirable path, through which the spectral optimisation should progress. 
The various spectral contours are chosen so, that the both the EVM and the 
PSD improves progressively as the spectral optimisation is progressing and 
local minimums are avoided* At this stage it is too early to specify these 
contours, but it is envisaged that after extensive characterisation an 
optimum contour set can be built. 

The stability and the convergence speed of the linearisation loop are 
critical parameters. Hence, it is crucial that the noise is suppressed as far as 
possible. Furthermore, error in the PSD relationship may contain several 
local minima. Random search must not be used to jump through these 
locations because it would lead to gross mask violations. 

The spacing of the intermediate masks should be close enough to produce 
a sharp error function between the measured PSD and the current contour 
thus steering towards the global minima. Initially the PSD is compared 
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with the first (i.e. the worst spectra) contour and the error is minimised 
within the iterative lineariser. As soon as the error becomes zero, the 
iteration continues with the second contour within the set. The process 
eventually reaches the target contour, which is then continuously held. One 
possible method for building intermediate spectra contours is illustrated by 
the following flowchart of figure 9. 
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CLAIMS 

1 . A method of linearising the output signal of a component having an 
input signal comprising; 

a) taking an initial estimated representation of the component 
inverse transfer characteristic; 

b) manipulating said representation; 

c) pre-distorting said input signal by applying values of said 
representation to the input signal; 

d) determining the error in terms of linearity between the input 
signal and the output signal; and 

e) repeating steps b) to d) such that the error is minimised. 

2. A method as claimed in claim L wherein said representation is in the 
form of a series of points; and wherein step b) comprises moving one or 
more of the points at each iteration. 

3. A method as claimed in claim 2 wherein step b) comprises moving 
one point at a time in a sequential order of increasing or decreasing 
scale. 

4. A method as claimed in claim 1 wherein said curve is represented as 
a series of basis functions in a NURBS curve representation; and 
said manipulation comprises scaling of the individual basis 
functions. 
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5. A method as claimed in any preceding claim wherein step d) 
comprises determining the error between the measured power 
spectral density and the target spectral density. 

6. A method as claimed in claim 5 wherein step d) minimises the area 
between the target power spectral density and the measured power 
spectral density across the full range of the spectrum 

7. A method as claimed in claim 5 wherein step d) minimised the error 
between target power spectral density and the measured power 
spectral density within a specific region of the power spectrum. 

8. A method as claimed in any preceding claim wherein step a) 
assumes no pre~distortion. 

9. A method as claimed in any preceding claim wherein said 
component is an amplifier. 

10. A system to linearise the output signal of a component having 
comprising: 

a) means to take an initial estimated representation of a component 
inverse transfer characteristic; 
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b) means to manipulate said representation; 

c) means to pre-distort said input signal by applying values of said 
representation to the input signal; 

d) means to determine the error in terms of linearity between the 
input signal and the output signal; and 

e) means to repeat steps b) to d) such that the error is minimised. 

11. A system as claimed in claim 10 wherein said representation is in 
the form of a series of points; and wherein said means to manipulate 
comprises means to move one or more of the points at each 
iteration. 

12. A system as claimed in claim 1 1 wherein manipulation means 
comprises means to move one point at a time in a sequential order of 
increasing or decreasing scale. 

13. A system as claimed in claim 10 wherein said curve is represented 
as a series of basis functions in a NURBS curve representation; and 
said manipulation means comprises means to scale the individual 
basis functions. 

14. A system as claimed in any of claims 10 to 13 wherein said means 
to determine error comprises means to determine the error between 
the measured power spectral density and the target spectral density. 
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A system as claimed in any of claims 10 to 14 wherein said 
component is an amplifier. 
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